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Abstract 

Twisted photons (i.e. photons carrying non-zero orbital angular momentum) are 
well-known in optics. Recently, it was suggested to use Compton backscattering off an 
ultra-relativistic electron beam to boost optical twisted photons into the high energy 
range. However, only the case of strictly forward/backward scattering has been studied 
so far. Here, we consider generic kinematic features of processes in which a twisted 
particle scatters with non-zero transverse momentum transfer. 

1 Introduction 

In perturbative quantum field theory we assume that interaction among the fields can be 
treated as a perturbation of the free field theory. This perturbation leads to scattering be- 
tween asymptotically free multiparticle states, which are usually constructed from the plane 
wave one-particle states. This choice greatly simplifies the calculations and represents a very 
accurate approximation to the real experimental situation in virtually all circumstances. How- 
ever, one can, in principle, choose any complete basis for the one-particle states other than the 
plane wave basis, provided that it is still made up of solutions of the free field equations. Such 
states can carry new quantum numbers absent in the plane wave choice and, if experimentally 
realized, they can offer new opportunities in high-energy physics. 

Thanks to the recent experimental progress in optics, it is now possible to create laser 
beams carrying non-zero orbital angular momentum (OAM) [I], for a recent review see [2]. 
The lightfield in such beams is made of states which are non-plane wave solutions of the 
Maxwell equations. Each photon in this lightfield, which we call a twisted photon, carries a 
non-zero OAM quantized in units of h. Several sets of solutions have been investigated, such as 
Bessel beams or Gauss-Laguerre beams, but in all cases the spacial distribution of the lightfield 
is necessarily non-homogeneous in the sense that the equal phase fronts are not planes but 
helices. Such states form a complete basis which can be used to describe the initial and final 
asymptotically free states. Moreover, it is the basis of choice for experimental situations when 
the initial states are prepared in a state of (more or less) definite OAM. 

Twisted photons have been produced in various wavelength domains, from radiowave [3] 
to optical, with prospects to create a brilliant X-ray beam of twisted light in the keV range 
jl]. Very recently it was suggested to use the Compton backscattering of twisted optical 



photons off an ultra-relativistic electron beam to create a beam of high-energy photons with 
non-zero OAM [51 [6]. The technology of Compton backscattering is well established, and 
the high-energy electron beams and the OAM optical laser beams are already available. The 
suggestion of [HI E] paves the way for the twisted photons, and the twisted states in general, 
into the high-energy physics. The wealth of new physics opportunities related to this new 
degree of freedom is yet to be understood (see [7] for initial steps in this direction). 

In this paper we begin this exploration by focusing on a technical question of how to 
calculate the scattering of a one-particle twisted state off plane wave states: 



Here systems X and X' are one- or multi-particle states described by plane waves, p and p' 
being their respective total momenta. We consider the case of non-forward (and non-backward) 
scattering, that is, when the transverse momenta of systems X and X' with respect to the 
axis used to define the twisted states are different, p ^ p'. 

In this work we focus on the general kinematical properties of the scattering matrix which 
hold for all processes of type ([1]). The purpose of this exercise is twofold. The main goal is to 
develop the formalism of treating the scattering processes with twisted particles in the initial 
and final states. The secondary goal is to extend the original calculation of [51 E] of the strictly 
backward Compton scattering of twisted optical photons to non-zero transverse momentum 
transfer. The basic questions are: how to write the cross section of this process, and how the 
parameters of the final twisted photon depend on the momentum transfer. 

Since the kinematical details turn out to be rather unconventional, we describe them in 
a pedagogical manner with the simplest possible process: the decay of a massive twisted 
scalar into two massless scalar. In this way we separate the effect of non-homogeneous spatial 
distribution from the possible effect of non-trivial polarization states, which the twisted photon 
can have and which we postpone for a future study. We stress however that the kinematical 
features we explore with this example are pertinent to all the scattering processes of type ([TJ). 
Additional process-specific properties arise on top of these kinematical features. 

The paper is organized as follows. In Section [5] we give an introduction into twisted particles 
and describe some of its properties. In Section E] we describe the general problem we tackle: 
scattering of a twisted particle accompanied by a non-zero momentum transfer. We pinpoint 
the key quantity to analyze, and then, in Section HJ we investigate this quantity in detail with 
a simple pedagogical example. Section |5] contain discussion of the results obtained and in the 
final Section we draw our conclusions. In Appendix we derive some properties of the Bessel 
functions used in the main text. 

2 Describing twisted states 
2.1 Spatial distribution 

As mentioned in the introduction, we focus in this paper on twisted scalar particles only. In 
this Section we follow essentially [3 E] • 

We represent a state with non-zero OAM with a Bessel beam-type twisted state. This 
is a solution of the wave equation in the cylindric coordinates with definite energy u and 
longitudinal momentum k z along a fixed axis z, definite modulus of the transverse momentum 
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|k| (all transverse momenta will be written in bold) and a definite ^-projection of OAM. If 
the plane wave state |PW(k)) is 

\PW{k)) = e - wt+ikzZ -e ikr , (2) 
then a twisted scalar state | k, m) is defined as the following superposition of plane waves: 

\ K: m) = e - Mfc ^|^a Km (k)e 4kr , (3) 

where a Km (k) = (-t) m e im ^V2^ 6 ^~ ^ . (4) 



In the coordinate space, 

pim<f> r 

\K,m) = e - Mk * z ■ Vw(r) , xp Km (r) = —=y/iiJ m (ivr) . (5) 

V Z7T 

Here, following [5] we call k the conical momentum spread, m is the ^-projection of OAM, 
and the dispersion relation is k^k^ = u 2 — k 2 z — n 2 = M 2 . We note in passing that the average 
values of the four-momentum carried by a twisted state is 

(h») = (co, 0, k z ) , (6) 

so that (k' 1 ){k l j) = M 2 + k 2 , which is larger than the true mass of the particle squared. 
The transverse spatial distribution is normalized according to 

A^M*»M = /nfrJ.WJUW = - «') ■ (T) 

The plane wave can be recovered from the twisted states as follows: 

|JW(k=0)) = limJ^I/c.O), (8) 

K->0 V ft 

/2tF + °° 

|PW(k)) = J— i m e~ im<l,k \K,m} , re=|kj_|. (9) 

m=— oo 

If needed, these two cases can be written as a single expression: 

\PW(k)) = lim W— i m e~ im<pk \K,m) . (10) 

^ m=— oo 

From these expressions one sees that the twisted states with different m and k represent 
nothing but another basis for the transverse wave functions. 

2.2 Density of states 

When calculating cross sections and decay rates, we need to integrate the transition proba- 
bility over the phase space of the final particles. When calculating the density of states, we 
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consider large but finite volume and count how many mutually orthogonal states with pre- 
scribed boundary conditions can be squeezed inside. In the present case due to the cylindrical 
symmetry of the problem, we choose a cylinder of large radius R and length L z . In the case 
of plane waves we have 

dn PW = tiR L z 3 . (11) 

The full number of states with transverse momenta up to and longitudinal momenta \k z \ < 
k z0 is k z0 L z ■ R 2 nil 'Air. 

To count the number of twisted states \k, m) in the same volume, we specify the boundary 
condition, e.g. ip Km (r = R) = 0, which makes k discrete such that KiR is the i-th root of the 
Bessel function J m . Then we note that the position of the first root of the Bessel function 
J m (x) is always at x > m, and as m grows x — > m. For a given k, the maximal m for which 
the wave can still be contained inside the cylindrical volume is m max = kR, which has a very 
natural quasiclassical interpretation. 

If m is small and not growing with R, then one can use the well-known asymptotic form 
of the Bessel functions to count the number of states: 



Rdn L z dk z Am 



dn tw = . (12) 



Here, Am is written instead of just 1 to signal the presence of a discrete running parameter 
m. 

If m is not restricted to small values, this asymptotic form of J m (x) cannot be used since 
it requires m 2 < x. Instead, the so-called approximation by tangents can be used, which gives 
the following density of states: 

-dn Am L z dk z 



dn tw = Vrn 2 max -m 2 - — ^ . (13) 



In the limit m m max = kR this expression reproduced ( fT2j) . Alternatively, one can calculate 
the radial part of the density of states via the adiabatic invariant as suggested in [6]. The 
number of radial excitations n r for a fixed m is 



f R k r {r)dr j m 2 

n r = — — , k r (r) = dK 2 - — . (14) 

The density of states is then given by 

dn r = -^dn = a/ m 2 — m 2 — . (15) 

dK " KTT 

One important remark is in order. Effectively, switching from the plane wave to twisted 
state basis for the final particles implies replacement 



m 2 , Am , . 

d 2 k^4./l — ndn . (16) 



m max mr 



Note that the contribution of each "partial wave" with a fixed m vanishes in the infinite 
volume limit as 1/R. However, the number of partial waves grows oc R, and in order to get 
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a non-vanishing result for a physical observable, one must integrate over the full available m 
interval up to m max . This holds even if the transverse momenta stay small, and is related to 
the fact that the plane wave contains contributions from all impact parameters with respect 
to any axis non-collinear to its momentum. 

Another expression one needs for the probability calculations is the normalization constants 
for the one-particle states. A usual plane wave one-particle state is normalized to 2E ■ V; to 
renormalize it to one particle per the entire volume, the plane wave should be multiplied by 
N PWl with 

N PW = ^y> V = rrR"L z . (17) 
For a twisted states the corresponding normalization factor N tw is 



1 TTK 



2E V m max ~ m2L z ' 

which in the small-m case simplifies to 

tw 2E RL Z ' 1 } 

which was also derived in [6]. Note however that even in the general case the product of 
the normalization constant squared and the density of states for each final twisted particle is 
simplified as 

„ t9 , d.Kdk z Am , , 

= inrsr- (20) 

3 Non-forward scattering of a twisted state: generic fea- 
tures 



k, m k', m' 




Figure 1: A generic process of non- forward scattering of a twisted state off plane waves. 

Consider again the generic process (JT|) which is schematically depicted in Fig. [TJ Here we 
have one twisted particle in the initial and one in the final states. Since we deal here with the 
kinematical features of the process, it is inessential whether these particles are of the same type 
or not. These two particles are described by the states \k, m) and \k', w!) defined with respect 
to a common z-axis. All the other particles are assumed to be plane waves, which allows us 
to define the total transverse momenta p and p' of the systems X and X', respectively. The 
transverse momentum transfer, q = p — p' is also well-defined, and we can consider two cases: 
(i) strictly forward scattering, q = 0, (2) non-forward scattering, q ^ 0. Note that we use the 
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term "forward" for any process with zero transverse momentum transfer, i.e. both for strictly 
forward and strictly backward scattering. 

Suppose that we know the .S-matrix for the same process with all particles represented by 
the plane waves. If the would-be twisted particles in the initial and final states have momenta 
k and k', the plane wave S"-matrix has the familiar form: 

S PW = i(27r)W 4 \k + p-k'-p')- M(k 2 , k' 2 , <f> w ). (21) 

The invariant amplitude Ai depends among other on the transverse momenta squared k 2 and 
k' 2 as well as on the azimuthal angle <pkk' between them, which is explicitly indicated in (121 p . 
Any other azimuthal angle, say, between k and p, can be written as <fik p = 4>k q + 0g P , and 
therefore it is expressible via (fikk> and internal angles in the systems X and X' . 

To get the ^-matrix for the twisted particle scattering process (jTJ, one just needs to apply 
fll]) to the initial and final states [5]: 

f d 2 k d 2 k' 

Stw = J j^^j^T^a* K , m ,(k')a Krn (k) S PW ■ (22) 



The delta-functions present in (I22j) fix k 2 = and, as we will see later, they specify 

the angle <fikk' up to its sign. If the plane wave matrix element is not sensitive to the sign of 
the angle <pkk'i Ai{n 2 , k' 2 , <fikk') — A4(k 2 , k' 2 , —4>kk'), then it can be taken out of the integral, 
and we are left with the following transverse process-independent master integral: 

/d 2 k d 2 k' 
(2ir) 2 (27r)2 a ^ m, ( k ') aKm ( k ) ( k ~ k ' + ^ • ( 23 ) 

If the matrix element depends on the sign of <pkk', then it can be decomposed into a symmetric 
and antisymmetric terms 

M = M{<f) k k') + M(-<p k k>) M = M(<f> kk ') - M{-<j> kk >) ^ 

then in the twisted scattering amplitude Ai s will be accompanied by the master integral 
fl23|) . while Ai a will enter together with a slightly modified version of the master integral, 
i m ,m'(K, K,',q), in which the integrand of ( l23l) is multiplied by the sign function e((f>kk')- In 
fact, two interfering plane wave amplitudes appearing here might lead to novel observable 
effects in twisted particle cross sections, see details in [7]. 

In the next Section we compute the master integral ( 123]) and explore the singularities of |X| 2 
which appear in the cross section/decay rate calculations. However, let us briefly summarize 
the results right away: 

• In the strictly forward case, q = 0, X mm i oc 5(k — K')S mtm >, that is, the quantum numbers 
of the twisted state are transferred to the final particles without any change. 

• In the non-forward case, 0, a distribution over k' arises, and k' can be large if the 
momentum transfer is large. 

• At any non-zero momentum transfer instead of <5 mjm / we observe a distribution of m' 
over the entire possible range of values, —m' max <m'< m' max , where m' max = k'R — > oo. 

• \I\ 2 contains singularities at the border of the kinematically allowed region, which must 
be carefully dealt with. 
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4 Decay of a twisted scalar 



We would like to explore the integral ( 123|) in the context of the simplest possible problem: 
decay of a twisted scalar particle with mass M into a pair of massless distinguishable particles 
due to the cubic interaction g ■ <3>0i</> 2 . To make the presentation more pedagogical, we will 
first calculate the decay rate when both particles in the final state are plane waves, then for 
the plane wave plus twisted final state, and finally for the case when both final particles are 
twisted. We will calculate the decay width in the center of mass frame defined by k z = 0. This 
is not the true rest frame because due to the transverse motion a twisted particle is never at 
rest. 



4.1 Two plane waves 

For future reference, let us first recall the calculation for the standard case when all the 
particles including the initial one are plane waves. The S'-matrix is given as usual by S = 
i{27r) 4 S^(p-k 1 -k 2 )-g. When squaring the delta-function, we use the standard prescription 

VT 

[6®(p -h- k 2 )] 2 — ► 6^(p -h- k 2 ) ■ 5^(0) = 5^{p -h- h)-^— . (25) 

(2tt) 4 

We also use the plane wave normalization for all the particles, The decay probability per unit 
time for a particle at rest is: 

(27T) A g 2 S^(p-k l -k 2 )VT /Ar2 . n . . 

dT = V. [N 2 PW f ■ dn PW {k x )dn PW {k 2 ) 



g 2 5(M-ujx-u) 2 ) 3 



(2vr) 2 

so that the total width is 



d 6 ki , (26) 



<f (27) 



16vrM 

Now we repeat this calculation for the initial twisted state \K,m), while keeping the plane 
wave basis the final particles. The S'-matrix is 

S = i{2n fg5{E-uj l -uj 2 )5{k lz + k 2z ) J ^a Km (k)^(k - k x - k 2 ) (28) 

= l (27:Yg5{E-u 1 -u 2 )5{k lz + k 2z )^^e im ^ 1 ^_ 12) , (29) 

where k\ 2 = ykf + k| + 2|ki||k 2 | cos(0i — (f) 2 ) and 12 is the angle of the 2D vector k i2 w.r.t. 
some axis x. The phase factor is inessential and m disappears in the decay rate. 

When squaring the above expression, we encounter the square of 5(k— ki 2 ), which is treated 
as in [3 |6] : 

[5(k - k 12 )f = 5{k - k 12 )5(0) -> 6(k - k 12 )- . (30) 

7T 

This prescription comes from the observation that at large but finite R and at k = k\ 2 the 
divergent integral in (128]) is regularized as 

poo pR p 

6(0) = / rdr[J m {Kr)] 2 / rdr[J m {nr)} 2 w - . (31) 

JO JO 7T 
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Thus, with the normalization factors plugged in, the decay rate has the form 



dT 



(2tt)V 



-5(E -ux- u 2 )5(k lz + k 2z )TL, 



5(K-k 12 )R 1 7T Vd?k x Vd?k 2 



8Euiu 2 ■ T 
g 2 8(k - k 12 ) 5(E -ux- u 2 ) 



7T V 2 RL Z (2tt) 3 (2tt) 



(27T) 



/V 



8EuxU) 2 



dk z d 2 ki d 2 k 2 . 



For the transverse integral we write 
5(K-k 12 ) 



K 



2 j d(j) 2 5 [k 2 -k\-k 2 2 - 2\kx\\k 2 \ cos(0i-0 2 )] = i 



where 



A = - W2(kfk| + K 2 k 2 + H?k 2 2 ) - - k\ - k\ 



(32) 



(33) 



(34) 



is the area of the triangle with sides k, |kx| and |k 2 |. The angular integral can be non-zero 




Figure 2: The allowed kinematical region of the values of |ki| and | k2 1 for a fixed k defined 
by the "triangle rules" fl35|) . 



only if such a triangle can be formed, that is if |ki| and |k 2 | satisfy the "triangle rules": 

« < |ki| + |k 2 | , |k x | < k + |k 2 | , |k 2 |</c+|ki|. (35) 

The allowed values of | ki | and | k 2 1 for a fixed k are shown in Fig. [21 Note that the angular 
integral (1331) receives equal contributions from two points: 



} i — 02 = ±^12 ; where 5 12 = arccos 



k 2 - k 2 - k 2 



As usual, the energy delta function can be killed by the k z integration 

S(E-ux-u 2 ) 2 



dk. 



UxU-2 



Ek* 



where 



2E 



E A + kf + k 4 - 2(E 2 k 2 + E 2 k\ + k 2 k 2 



(36) 

(37) 
(38) 



The decay rate becomes 

g 2 1 NjkjJ^MN ^ 
dL ~4^4E*k* z A • [M) 

The integration region over |ki| and | lea | is defined by the requirement that, in addition to 
fl35|) . the longitudinal momentum k* z is well-defined, which cuts the rectangular shape shown 
in Fig. [2] from above. It is conveniently described with variables 

|ki| — |k 2 | |ki| + |k 2 | _ E 

x = , z = , x E [-1, 1J , z G [1, z max \ , z max = — > 1 . (40) 

Kj Kj K 

In these variables 



A = ^vV-l)(l-* 2 ), K = jJ(l~iP)(l-^), (41) 



max / \ max 



and the decay rate takes form 

dT = -4 ^ ~ x2)dzdX m . (42) 

[(^-^)(^L,-^)(^-i)(i-^ 2 )] 1/2 

This integral can be taken exactly, and it gives the decay width of the twisted scalar particle 

9 2 (43) 



IQtxE 

which is a very natural result. In the limit k — > 0, we recover the plane wave decay width 

urn 



4.2 Twisted state plus plane wave 

Let us now describe the final state as twisted state plus a plane wave: 

|«, m) ->• |ki, mi) + |/W(k 2 )) . (44) 

Since the full decay width cannot depend on the basis we choose for the final particles, we 
must recover the same result ( )43l) in this basis. In addition to that, we also want to know how 
the final twisted state parameters Ki and mi are related to the initial state parameters k and 
m. 

4.2.1 Forward case 

The S- matrix is now expressed in terms of the master integral ( 1231) : 

S = i(2ir) 4 g5(E u 2 )5(k zl + k z2 ) ■ X m>mi (K, k 2 ) . (45) 

Let us first calculate this integral in the strictly forward case k 2 = 0: 

X m , mi (^ 1; 0) = ^ / d 2 kd 2 kie^-^^ ^^^4^ ^)(k-k x ) 



^— y ##ie^ m ^ -25(k 2 -k?)5(0-</»i) 

— ^5(K - Kl)5 m ,mi , (46) 



which was first obtained in [5] . This result implies that the twisted state quantum number are 
transferred from the initial to the final twisted particle without any change. The differential 
decay rate is 

, r 9 2 6{E - uj x - u 2 )S{k lz + k 2z ) \ 7 ji ,3, 

(2ti) z 8EU1CU2 

- 3 ' « B -*»-«>^-*gL, (47) 



(2tt) 2 8£wiw 2 (27r) 2 4E 2 fc 
and we stress that this result is applicable only at k 2 = 0. 

4.2.2 Non-forward case: the master integral 

Now we consider the non-forward case. The master integral can be rewritten as 

X m>mi (K, Kl ,-k 2 ) = J^ys^ J d^e^™^' -^(k-ki-ka). (48) 
There are two ways to look at this integral. First, we can rewrite 

<5 (2) (k - ki - k 2 ) = J dh e trk ~ trkl ~ irk2 (49) 

and represent the integral as 

2m,mi(«,«i,k 2 ) = y/KKi- J rdr J to (kt) J mi («ir) J m _ mi (|k 2 |r) , (50) 

where 2 is the azimuthal angle of k 2 . Note that although we are considering the case with 
two twisted particles (one in the initial and one in the final state), an integral over three Bessel 
function arises automatically. 

The master integral can be also calculated by a direct integration over angles in fj48|) . We 
rewrite the delta-function as 

5 2 (k-k!-k 2 ) = 25[K 2 -(k 1 + k 2 ) 2 ]5(0-0k 1+k2 ) 

= 25{k 2 - k\ - k 2 , - 2/ci|k 2 | cos(0i - 2 )] <5(0 - 0k 1 +k 2 ) • (51) 

One sees that the integral can be non-zero only if it is possible to form a triangle with sides 
k, K\, and |k 2 |; thus, K\ and |k 2 | satisfy the same triangle rules as in (135]) . Fixing k\ and |k 2 | 
means that the triangle can be formed only at two choices of relative azimuthal angles of k, 
k x and k 2 . Let us introduce 

/ /t 2 + /t? ko \ _ / /t 2 + ko /Kn\ 

0\ = arccos , o 2 = arccos — : , (52) 

\ 2kki J \ 2ft|k 2 | / 

Then, the two configurations of transverse momenta correspond to 

4> - 0i = ±£l , 4> - 4>2 = =F^2 , (53) 
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so that the signs of — <j>\ and — 2 are always opposite. The master integral then becomes 

= g^^'— ) ^° S|m ^' - . (54) 

where A is the same as in (134"1) with k 2 replaced by k 2 . Comparison of (1501) with (154")1 gives the 
result for the integral of the triple Bessel function product. See also [8] for some mathematics 
involved in evaluation of this and similar integrals. 

Let us also check the k 2 — > limit of the result (I54p . When |k 2 | <k, the distribution over 
K\ spans from k — | k2 1 to k + | 1c2 | , see Fig. [2j The angle 5\ — > 0, while <5 2 can still be arbitrary. 
However, in the k 2 — > limit only m — mi = term survives due to J m -m 1 (|k 2 |r), so that the 
cosine in (l5~4j) approaches unity. The analysis of A shows that 

lim = 2tt5(k - k x ) , (55) 

|ka|->0 A v ' 



so that one indeed recovers the strictly forward result for the master integral given in (14611 . 
Note that although we wrote simply |k 2 | — > 0, the exact definition of this limit is involved and 
is discussed in Section [5721 

Note also that the modified master integral X m m /(K, k', q) mentioned in the previous Sec- 
tion has the same form as (l54"j) but contains sine instead of cosine of rriiSi + (m — mi)5 2 . This 
modified integral does not enter our particular calculation since the amplitude is symmetric 
under the 4>kk' sign flip. 

4.2.3 Non-forward case: squaring the amplitude 

After integrating over k z , the decay rate can be written as 

dT = • |X m , mi («, «4, k 2 )| 2 ■ ^ d 2 k 2 . (56) 

Note that this decay rate is differential not only in k,x and k 2 but also in the discrete variable 
mi, the full decay width includes integrals over momenta and a summation over all possible 
mi's. 

A close inspection shows that the immediate integration over /«i or k 2 cannot be done 
due to singularities of |X| 2 along the boundaries of the kinematically allowed region shown in 
Fig. [2j In contrast to the plane wave case, the denominator now contains A 2 instead of just 
A. Therefore, in terms of variables x and z one encounters singularities of the form 

dx , f Zmax dz 

and 



-1 



1 — x 2 /1 z 2 



Clearly, this is an artefact of the infinite radial integration range. If instead we take R to be 
large but finite, we expect that a trick similar to (1301) should be at work, namely that after 
regularization |X| 2 would yield R times a less singular function. 

This trick does not seem to work for each mi separately. However, as we prove in Appendix, 
it works for |X| 2 summed over all possible mi. In the limit R — > 00 we obtain 



E 



mi=— 00 



.-, 1 Rk\ 1 



iX^O^k^l (57) 
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with the same A as before. The regularization parameter R then disappears from the result, 
and the decay rate reads 

rp _ 9 2 1 M^i |k 2 |d|k 2 | , , 

dT ~^A^K A • (58) 

Comparing ( 1581) with the previous results (H71) and ( 1391) leads us to two conclusions. 

• The transition from the strictly forward to the non- forward cross section/decay rate 
consists in replacement 

_1 (59 ) 

• The result ( 15 8p for a twisted particle in the final state coincides with the result (1391) for 
the case when both final particles are plane waves. 

Although these conclusions were drawn for the specific process we consider, the way it is 
derived suggests that this might be a universal feature for many (or all) processes of type (JTJ). 

4.3 Two twisted states 

For completeness, let us also recalculate the decay rate in the basis when both final particles 
are described by twisted states: and \n 2 ,m 2 ). We remind that all twisted states are 

defined with respect to the same common z axis. It turns out that this calculation closely 
follows the case of twisted state plus plane wave just considered. This is not surprising because 
the appearance of the triple Bessel integral highlights the fact that when two particles (one 
in the initial and one in the final state) are twisted, the third one is automatically projected 
from the plane wave onto an appropriately defined twisted state as well. 
The S"-matrix takes the form 

S = ig(27r) 3/2 5{E — ui — u 2 )5(k z i + k z2 ) y/K,K 1 n 2 J rdr J m (Kr)J mi (Kir)J m2 (K 2 r) , 

(60) 

and we again encounter the triple Bessel-function integral. The decay rate is written as 

dV= A^A&K- A • (61) 

This expression is identical to fl58|) up to the obvious replacement |k 2 | — >■ « 2 ; its integration 
over all conical momenta spreads Ki, k 2 gives again 



5 Discussion 

5.1 Distribution over Kj 

Let us first discuss what typical values of K\ and mi essentially contribute to the decay rate. 
The differential decay rate (1581) shows that at large |k 2 | 3> k the conical momentum spread 
K\ is limited to the interval from |k 2 | — k to |k 2 | + k with an inverse square root singularity at 
the endpoints. This singularity is integrable, so that the entire interval more or less equally 
contributes to the integral. Since |k 2 | can be as high as E, the total decay width is therefore 
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dominated by large K\ ^> k. In a more complicated process, the decay rate or the cross section 
will include the amplitude squared which can serve as a cut-off function. For example, in the 
Compton scattering one expects that K\ up to ~ m e will contribute to the cross section. 

5.2 Distribution over m\ 

Although we found a result for the decay rate summed over all mi, we can trace the main 
mi-region from the intermediate formulas. The result is that essentially all ni\ from minus 
to plus infinity are important for the decay rate, which is in a strong contrast to the strictly 
forward result m — mi. 

Indeed, the distribution over final mi can be seen in (1561) . where the master integral X m ^ m i 
is given by fl54l) . For generic transverse momenta, growth of mi leads to oscillations of the 
cosine function with constant amplitude. Thus, when averaging over the entire Ki interval, 
one can approximate cosine squared by 1/2, and the dependence on mi drops off. This result 
holds for any non-zero transverse momentum transfer | ki2 1 — > 0. 

Since the forward and non-forward mi-distributions are so dramatically different, a natural 
question arises whether there is a continuous transition from the non-forward to the forward 
scattering. The answer to this question involves an accurate treatment of two limits: [ k2 1 — > 
and R — > oo. Let us keep R large but finite, and set | k2 1 — > 0; then the transition is 
smooth. Looking at the triple-Bessel representation of the master integral f )50p with the upper 
limit replaced by R, one sees that the result will begin to significantly decrease only when the 
position of the first node of the last Bessel function J m - mi (|k 2 |r) falls outside of the integration 
range, that is for |k 2 | < |m — rrii\/R, where m-i ^ m. If m x = m, then at |k 2 | <C 1/R the last 
Bessel function can be approximated by the unity. Therefore, the limit |k 2 | — > taken in f )55|) 
implies that 

I k 2 1 — > and R — > 00 provided that | k 2 1 R -C 1 . 

If instead R — > 00 at fixed |k 2 |, then the transition of non- forward to forward results is 
discontinuous at |k 2 | = 0. 

In [5] it is claimed with the specific example of the Compton cross-section that if the 
transverse momentum transfer is small compared to k (in our notation, finite |k 2 | <C k at 
infinite R), then the mi-dependence has a narrow distribution peaked at mi = m. Our 
analysis shows that this cannot be true, and this conclusion of [5] is likely to be the result of 
an incorrect small- |k 2 | approximation for X mjmi (K, m, k 2 ). 

5.3 Orbital angular momentum vs. orbital helicity 

Looking back at the formalism used, we can conclude that our result that all mi essentially 
contribute to the decay rate/cross section just reflects the unfortunate choice of the same 
common axis z for all the twisted states appearing in the process. It does not give a clue of 
how twisted the final particles are with respect to their own propagation axes defined by their 
average values of the 3-momentum operator. Indeed, even a simple non-forward plane wave 
when expanded in the basis of twisted states contains all partial waves, see Nevertheless 
it carries a zero orbital angular momentum with respect to its own direction of propagation. 
Therefore, it appears that a more physically reasonable quantity is the "orbital helicity", 
projection of orbital angular momentum on the axis of motion. The relevant question is then 
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how this "orbital helicity" , not the OAM with respect to a fixed axis, is transferred from the 
initial to the final twisted state. We postpone this question for future studies. 



6 Conclusions 

Photons carrying non-zero orbital angular momentum (twisted photons) are well known in 
optics. Thanks to the recent suggestion [51 [6] to use the Compton backscattering of optical 
twisted photons off a high-energy electron beam, twisted photons are now entering the high- 
energy physics. If the new degree of freedom they offer is effectively realized in experiment, 
twisted photons can bring novel opportunities to particle physics. 

In this paper we took a first look at the non-forward scattering of a twisted particle off a 
system of plane waves: 

twisted + X(p) twisted' + X'(p'), p7p, 

and investigated universal kinematical features pertinent to such processes. 

With the simple example of the decay of a scalar twisted particle, we evaluated the master 
integral appearing in such processes, discussed its singularities and explained how to write 
the cross section/decay rate. These results can be now used, for example, to investigate the 
Compton backscattering of the twisted photons in the non-forward region, which was missing 
in the original suggestion [5j E] . 

Discussing the results obtained, we came to the conclusion that a more physically motivated 
quantity to describe a twisted state would be the orbital angular momentum projection not 
on the reaction axis but on the direction of the outgoing twisted particle. Incorporation of 
this "orbital helicity" into the present formalism remains to be done. 
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A Regularization of \I\ 



Here we calculate the large-i? behavior of the mi-sum of the squares of the triple-Bessel 
integral: 



E 



mi- 



rdrJ m (Kr)J mi (Kir)J m - mi U 2 r 



(62) 



which appears in the decay rate (1561) . Evaluation of the integral itself with R — > oo performed 
in the main text shows that it can be non-zero only if k, k\, k 2 satisfy the triangle rules (135 p . 
i.e. a triangle with these sides can be constructed. Since m describes the initial state, we take 
it small and not growing with R: m <^ m\ max = KiR, while mi can extend up to m\ max . The 
final Ki, k>2 can be much larger than k. 
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Since the expression (1621) is regularized with large but finite R, the summation and inte- 
gration can be interchanged: 



rdrr'dr' J m (Kr)J m (Kr') ^ 4, far)J m - mi far)J mi far')J m - mi (« 2 r / ) . (63) 



mi=— mimu 



Thanks to the properties of the Bessel functions, only mi's up to mm(K 12 r, /Ci^r') are ef- 
fectively contributing to this sum; for larger mi the Bessel functions strongly decrease. But 
r, r' < R, which means that the limits on the summation can in fact be safely extended to the 
infinity. Then, the sum of the product of four Bessel functions is treated in the following way: 



^ Jmi fa r ) Jm-mi far) 4i far') 4- mi far') 
oo 
+00 

+00 ^ „2tt 

V — / dae iimi ~ m '^ a J mi far)J m - mi far)J m >far')J m - m >far') 

+00 

]T e- iro i a J m /(K 1 r , )J m _ mi (K 2 r'; 

m\—— 00 

The first sum in the square brackets is calculates as follows: 



(64) 



m\=— 00 

+00 



mi ,771]=— 00 



— / da 

2tt 



+00 



e imia J mi far)J m _ mi far) 

,mi=— 00 



e imia J mi (Kir)J m _ mi (/€ 2 r) 

mi=— 00 

(-0" 



(2tt) s 
2tt 



>i a© 2 e 



IK\T COS <pl+lK2»" COS (/>2 



:} imia+imiipi+i(m—mi)ip2 



^ g«m0i giKir cos 4>l+ik2r cos{<f>\+a) 



The combination of angles and momenta inside the exponential can be expressed as 

«i cos 0i + k 2 cos(0i + a) = |ki + k 2 | a cos(0i + 50) , 
where 

k 2 sin a 



ki + k2| a = \ k\ + «| + 2/«i«2 cos a , tan< 



(65) 



(66) 



(67) 



K\ + k 2 cos a 

Geometrically, |ki + k 2 | Q is the norm of the sum of two vectors of moduli K\ and k 2 and the 
relative azimuthal angle a. Therefore, (165)) is 



(68) 



^ a Jm,far)J„far) = e M«-W j m (| kl + k2 | Q7 

mi=— oo 
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This expression can be viewed as a 2D generalization of the well-known addition formula for 
the Bessel functions 



+00 



m\=— 00 



Now, the second sum differs only by a — > —a (or, alternatively, complex conjugation) and 
r — > r'. Therefore, the summation (I64I) is simplified to 

1 f 27T 

— daJ^lkx + kaUr) J m (|k x + k 2 | Q r') . (69) 
2vr Jo 

Note that this integral involves only the Bessel functions of small order. We now plug this 
expression in (1631) and get 



y da J rdr J m (nr) J m (\k 1 + k 2 \ar) J r'dr' J m {nr') J m (|ki + k 2 | a r') (70) 



As usual, we extend the integration range in one of the integrals to infinity, which gives a 
delta-function, and then we use it on the second integral calculated up to R. The original 
expression fl62l) then becomes 

where A is, as always, the area of the triangle with sides k, k±, k 2 . 
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